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THE  ANALYSIS  FOR  THE  REGULATION  PERFORMANCE  OF  A  VARIABLE 
THRUST  ROCKET  ENGINE  CONTROL  SYSTEM 

by  Zhou  Xiwen 
Abstract 

This  paper  introduces  the  constituent,  operation  principle 
and  regulation  process  of  the  bipropellant  variable  thrust 
rocket  engine  for  a  variable  area  injector  that  uses  solenoid 
valve  control. 

According  to  the  constituent  and  practical  operation  process 
of  an  engine  control  system,  the  mathematical  models  of  each  link 
is  given  and  the  transfer  function  method  is  used  to  analyze  the 
dynamic  and  steady  performances  of  engine  control  systems.  The 
use  of  concrete  numerical  calculations  further  shows  the  major 
factors  influencing  the  dynamic  performance  and  steady  errors  for 
the  engine  control  system. 

I.  The  Constituent  and  Operation  Principle  of  Engine  Control 

Systems 

This  paper  will  introduce  the  use  of  a  control  circuit  which 
possesses  a  certain  operating  frequency  and  changeable  pulse 
width  to  control  two  high  performance  solenoid  valves.  The  use 
of  the  flow  quantity  change  of  liquids  to  change  the'  pressure  of 
the  injector's  hydraulic  pressure  chamber  Will  cause  the  up  and 
down  movements  of  the  injector's  needle  valve  to  change  the  area 
of  the  flow,  control  the  flow  quantity  of  the  propellant  and 
realize  thrust  changes.  The  control  system  is  mainly  composed  of 
a  control  circuit,  solenoid  valve  and  thrust  chamber  with  a  var¬ 
iable  area  injector.  Fig.  1  shows  the  relationship  between  them. 
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Fig.  1  Systems  chart  of  the  operation  principle  of  the  engine. 

Key:  (1)  Control  voltage;  (2)  Delay  multivibrator; 

(3)  Solenoid  valve  drive  circuit;  (4)  Operational 
amplifier;  (5)  Pulse  vibrator;  (6)  Delay  multivibrator; 

(7)  Solenoid  valve  drive  circuit;  (8)  Control  circuit 
section;  (9)  Filling  device;  (10)  Electric  explosion 
valve;  (11)  Entrance  solenoid  valve;  (12)  Rubber  ring; 

(13)  Pressure  sensor;  (14)  Combustion  chamber;  (15) 
Spring;  (16)  Exit  solenoid  valve;  (17)  Needle  valve. 

The  control  circuit  is  composed  of  an  operational  amplifier, 
pulse  vibrator,  delay  multivibrator  and  solenoid  valve  drive 
circuit  etc.  The  operational  amplifier  carries  out  superposition 
operations  of  the  control  circuit  signals  as  well  as  the  combus¬ 
tion  chamber's  feedback  signals  and  supplies  voltage  error 
signals.  The  pulse  vibrator  gives  voltage  pulse  signals  with 
certain  frequency  and  amplitude  and  controls  the  delay  multivi¬ 
brator.  The  delay  multivibrator  controls  the  pulse  operation  of 
changing  width  based  on  the  frequency  of  the  pulse  vibrator  and 
size  of  the  operational  amplifier's  output  voltage  error.  These 
pulse  signals  which  possess  certain  frequency  width  variations 
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are  amplified  by  the  solenoid  valve  drive  circuit  and  drive  the 
solenoid  valve  pulses. 


The  two  solenoid  valves  are  connected  to  the  injector's 
hydraulic  pressure  cavity.  When  the  error  signals  of  the  control 
circuit  are  positive  voltage  signals/  the  control  entrance  sole¬ 
noid  valve  operates,  the  pressure  of  the  hydraulic  pressure 
cavity  rises,  the  needle  valve  of  the  thrust  injector  lifts  up, 
the  nozzle  area  of  the  injector  enlarges,  the  flow  of  the  propel¬ 
lant  increases  and  the  thrust  rises.  On  the  other  hand,  when  the 
error  signals  of  the  control  circuit  are  negative  voltage  signals, 
the  control  exit  solenoid  valve  operates  and  eliminates  the 
liquid  in  the  hydraulic  pressure  cavity  causing  the  pressure  in 
the  hydraulic  pressure  cavity  to  decrease.  When  the  injector's 
needle  valve  is  under  the  effects  of  the  spring  force,  the  flow 
area  of  the  nozzle  decreases,  the  flow  of  the  propellant  de¬ 
creases  and  the  thrust  also  decreases. 

The  control  working  substance  of  the  injector's  needle  valve 
uses  a  propellant  which  is  passed  through  the  entrance  solenoid 
valve  into  the  hydraulic  pressure  cavity.  Because  the  flow  of 
the  control  fluid  is  very  small,  it  passes  through  the  exit 
solenoid  valve  and  is  directly  discharged  into  the  atmosphere. 

We  can  also  design  a  system  wherein  the  control  fluid  returns 
again  to  the  incendiary  agent's  main  pipeline. 


tig.  2 

3 


Fig.  2  Diagram  of  the  structure  of  the  solenoid  valve. 

Key:  (1)  Exit  connecting  nozzle;  (2)  Needle  valve  seat; 
(3)  Needle  valve?  (4)  Spring;  (5)  Coil;  (6)  Entrance 
connecting  nozzle. 


We  can  see  from  the  above  that  this  type  of  variable  thrust 
rocket  engine  system  uses  the  change  of  the  size  of  the  injector's 
flow  area  to  change  the  flow  of  the  propellant  and  attain  thrust 
regulation.  The  propellant  supply  system  of  the  engine  does  not 
carry  out  regulation.  The  special  features  of  its  plan  are  that 
the  control  system  is  simple,  operation  is  reliable  but  regula¬ 
tion  precision  is  not  sufficiently  high. 


Fig.  3  Block  diagram  of  the  control  system. 

Key:  (1)  Control  voltage;  (2)  Control  circuit? 

(3)  Entrance  solenoid  valve;  (4)  Injector;  (5)  Combustion 
chamber;  (6)  Exit  solenoid  valve?  (7)  Pressure  sensor. 


The  operating  parameters  of  the  control  circuit  are  as 
follows: 

Source  voltage  of  the  control  circuit  V=28  [V] 

Control  voltage  of  the  circuit  Vc=0  ~  5  [V] 

Operating  frequency  of  the  pulse  oscillating  circuit 
f =50  ~  100  [Hi 

Z 

Output  voltage  of  the  feedback  pressure  sensor 
Vf=0  ~  5  [V] 

Output  voltage  of  the  operational  amplifier 
VA=0.1  ^  1.1  ~  2.1  [VJ 
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II.  Mathematical  Equations  of  the  Main  Links  of  the  Control 
System 

The  solenoid  valve  of  the  engine  control  system  is  controlled 
by  the  control  circuit  and  carries  out  variable  width  pulse  oper¬ 
ations.  Therefore,  we  establish  mathematical  models  for  the 
control  circuit  and  solenoid  valve.  During  mathematical  process 
ing,  the  Laplace  transformation  of  the  non-continuous  system  will 
be  encountered,  that  is,  the  problen\  of  the  Z  transformation. 

However,  the  control  system's  input  signals  -  control  voltage 

V  and  output  signal  -  thrust  chamber  pressure  p  are  continuous 
c  c 

signals  and  only  the  intermediate  links  are  pulse  signals.  In 
order  to  simplify  the  overelaborate  mathematical  deductions,  we 
establish  mathematical  models  for  the  main  links  such  as  the 
circuit,  solenoid  valve,  injector,  combustion  chamber  etc.  based 
on  the  practical  operation  and  regulation  state  of  the  engine 
system  as  well  as  the  processing  of  the  continuous  system: 

1.  The  Control  Circuit  Equation  Basic  Hypotheses: 

1)  The  pulse  operating  frequency  of  the  control  circuit  is 

fixed.  The  size  of  the  circuit's  operating  frequency  influences 

the  dynamic  and  steady  state  performances  of  the  engine  system. 

Moreover,  the  size  of  the  frequency  is  also  influenced  by  the 

performance  of  the  solenoid  valve  itself.  Tests  show  that  the 

operating  frequency  of  the  engine  system  which  we  introduced  can 

be  selected  within  the  range  of  50  1,000  H  . 

z 

2)  The  pulse  width  of  the  control  circuit  forms  a  direct 
ratio  with  the  voltage  deviation  (the  difference  between  the  con¬ 
trol  voltage  and  feedback  voltage  -  the  deviation  signal) . 

(Tests  prove  that  the  pulse  width  and  voltage  deviation  basically 
change  within  a  direct  ratio) . 

In  order  to  raise  the  dynamic  and  steady  state  performances 
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of  the  engine  system,  the  pulse  width  should  be  regulated  as 
wide  as  possible.  The  pulse  frequency  of  the  control  circuit 
and  the  pulse  width  do  not  influence  eachother  and  they  can  be 
separately  regulated. 

3)  The  leading  edge  starting  time  of  the  solenoid  valve  is 
fixed.  Because  the  solenoid  valve  is  an  inductive  unit,  the 
electromagnetic  force  must  overcome  the  fictional  force,  spring 
force,  damping  force  and  hydraulic  pressure  etc.  of  the  needle 
valve.  The  electromagnetic  force  is  in  direct  ratio  with  the 
size  of  the  electric  current  and  the  change  of  the  electric  cur¬ 
rent  with  the  time  rises  according  to  the  exponential  curve  and 
thus  forms  the  leading  edge  when  the  solenoid  valve  opens.  The 
operating  electric  current  curve  of  the  solenoid  valve  is  shown 
in  Fig.  4.  Starting  leading  edge  *ZJe=3  5ms,  the  opening  time 
of  the  needle  valve  'C^=0.5ms  and  is  called  the  effective 

pulse  width  of  one  pulse  of  the  solenoid  valve. 


Fig.  4  Operating  electric  current  curve  of  the  solenoid  valve. 

Based  on  the  above  analysis,  we  can  give  the  circuit  - 
solenoid  valve  equation. 

k-st  <» 

In  the  formula,  K.  is  the  circuit  gain  [S/V] ,  E  =V  -V_  is 

the  difference  between  the  control  voltage  and  feedback  voltage 

[V] ,  when  V  =0  then  E  -V  and  7?  .  is  the  effective  pulse 
r  o  c  d 
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width  [S] . 


After  the  engine  is  regulated  to  the  operating  state  and  the 
pulse  frequency  is  selected,  the  pulse  width  is  regulated  to  its 
maximum  width  and  K,  is  also  fixed. 

D 

For  example,  f=50Hz,  Eq=Vc=5V  corresponds  to  the  maximum 

pulse  width  and  when  %  +  *tv=20ms  (  7?  =5ms,  '2T.  =15ms,  7^,-  is 

e  b_  ,  e  ~  b  f 

deleted),  then  Kfe=  '£,b/Eo=3xlO~J  [S/V]  . 

The  pulse  width  equation  of  the  solenoid  valve  is 

(2) 

In  the  unit  time  of  the  solenoid  valve  there  is  a  iffective 
operating  time  called  the  solenoid  valve  effective  pu  r  operat¬ 
ing  coefficient.  It  is  indicated  by  H  (t) .  Naturally, 

7(0-/6(0  (3) 

In  the  formula,  f  is  the  pulse  operating  frequency  [Hz ]  of  the 
solenoid  valve  and  H (t)  is  the  effective  pulse  operating 
coefficient  (*)(t)  <U. 

2.  The  Solenoid  Valve's  Flow  Equation 

Basic  hypotheses: 

1)  After  the  solenoid  valve  passes  leading  edge  time  % 
within  a  pulse  operating  time,  the  needle  valve  instantaneously 
opens. 

Due  to  the  fact  that  after  the  solenoid  valve  passes  leading 
edge  time  7^=3  “’^'Sms  £he  time  that  the  needle  valve  is  open  is 
only  7?f=0.5ms,  ^  therefore  "£f  can  be  deleted. 

2)  Tests  prove  that  when  the  closing  time  -  trailing  edge  of 
the  solenoid  valve's  needle  valve  is  smaller  than  0.5ms,  it  can 
be  deleted  and  it  can  be  considered  that  the  needle  valve  is 
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instantaneously  closed. 


Based  on  the  above  hypotheses,  when  considering  the  structural 
measurements  of  the  solenoid  valve,  we  can  arrange  the  flow  equa¬ 
tion  of  the  entrance  solenoid  valve  and  exit  solenoid  valve. 

The  flow  equation  of  the  entrance  solenoid  valve  is: 

Q&O-itAyj  &ZP'-P„( 03 

In  the  formula,  M  is  the  flow  coefficient  and  A  is  the  flow 
2 

area  [mm  ]  of  the  solenoid  valve  which  is  calculated  from  the  geo¬ 
metric  dimensions  of  the  solenoid  valve  (see  reference  t2]): 

A  =  2.226(D-  )/2)‘ 

(S'  is  the  maximum  interval  [mm]  of  the  solenoid  valve’s  open 
needle  valve: 

3=0.  l~Q.4mm  • 

D  is  the  diameter  [mm]  of  the  solenoid  valve's  exit,  P  is  the 

2  ° 
fluid  pressure  [kg/mm  ]  and  P„(t)  is  the  hydraulic  pressure 

2  H 

cavity's  pressure  [kg/mm  ]  of  the  thrust  chamber's  head. 

The  exit  solenoid  valve's  flow  equation  is: 

QH0-iiaJ22.Ph(.0 


Fig.  5  Geometric  measurements  of  the  solenoid  valve. 
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The  fluid  of  the  exit  solenoid  valve  directly  discharges 
into  the  atmosphere. 

The  solenoid  valve  is  pulse  operated  and  because  of  this  the 
practical  flow  of  the  solenoid  valve  should  introduce  the  effect¬ 
ive  pulse  operating  coefficient  0  (t) .  Therefore,  the  flow 
equations  of  the  entrance  solenoid  valve  and  exit  solenoid  valve 
are : 

Q.CO-’jCOMv^-CS.-SkCO] 

0,(0='j(0mV4“Ph(0’ 

We  substitute  the  (t)  and  b(t)  of  equations  (2)  and  (3) 
into  the  above  equation  and  let 

KQ-fK>nA,J&-  (O 

Then,  the  above  flow  equation  becomes: 

Q1(0-A'«£(0v'P.-P*(0 

■  QrCO-KoEitWPjj) 

We  know  from  the  operating  principle  of  the  engine  system 
mentioned  above  that  when  the  deviation  of  the  control  voltage 
and  feedfack  voltage  is  E{t)  ^  0,  the  control  entrance  solenoid 
valve  operates.  On  the  contrary,  when  E(t)  ^  0,  the  exit  sole¬ 
noid  valve  operates,  that  is,  when  the  two  solenoid  valves 
operate,  which  one  operates  depends  on  whether  the  deviation 
signal  E(t)  is  positive  or  negative. 

In  order  to  use  one  equation  in  the  control  circuit  to  indic¬ 
ate  the  solenoid  valve's  flow  equation,  it  is  necessary  to  employ 
two  functions.  Letting: 

,(e)a2MT[i£(,)i+£(0]  (5) 

l _ 

y(c)a*2|E(t)l  CISCO! -£«)] 


(6) 


Naturally, 


E(t)  >  0,  u (e) =1 ,  v (e) =0 ; 
E(t)<  0,  u  (e)  =0  ,  v  (e)  =1. 
and  then  the  solenoid  valves'  flow  equations  are: 


Q1(t)*KQu(e)E(t)  Aj Po-PH(tj 

(7) 

Q2(t)=K0v(e)E(t)  VPH(t) 

(8) 

These  are  binary  nonlinear  equations  and  based  on  the  design 
requirements  used  for  the  engine,  the  regulation  amplitude  of 
each  thrust  is  relatively  small  and  the  regulation  changes  of  the 
thrust  are  transitionally  smooth.  We  can  use  the  Taylor  series 
to  linearize  the  two  equations.  For  the  entrance  solenoid  valves 
it  is : 


0,(0 


.30,(0  |  £•(»)— dQ»(0  j  \Pf.(t ■) 

~mn  ,  dP„w\, 


In  the  formula: 


Within  this, 


Letting 


30,(0 
3£(  0 


~K0u(e)JP.-PH,  ~KqKI 


WO  =KQ  l-u^E.lJP.-P^Y'-K'K: 
dPH(,*)  .  2 


q-u(e)  -V  VPHt 

(9) 

Kl=  i  u(e)Et 1  Vvv1'1 

(10) 

Ki=KQKi 

(11) 

K" =K.K" 
u  Q  i 

(12) 

Q1(t)=K’E(t)-iT  A  PH(t) 

(13) 

then 

We  obtain  the  flow  equation  of  the  exit  solenoid  valve  in 
the  same  way: 


q,(0-k;e(0+k:ap,(0 
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In  the  formula: 


k:=kqk'. 

OS) 

k:=*kqk; 

(16) 

t—  y(e)>/  P  m 

(17) 

(e)E,CVP„t 

(18) 

Equations  (13)  and  (14)  are  indicated  by  one  equation: 


Q  (t)  -Qx  (t)  -Q2  (t)  -  (K^-IT )  E  (t)  -  (K"+lT )  A  PR  (t) 
Letting 


K  =K ' -K ' 
e  u  v 


kph=k;+k 


II 

V 


then 

Q  ( t )  =K  E  ( t )  -K  <dPH(t)  (19) 

e  PH  H 

Naturally,  when  E(t)  >0,  Ke=K^,  Kp  ,  Q (t) -Qx (t) ;  E(t)  <  0, 
Ke=Kv  '  KP  =Kv'  Q(t)=~Q2(t)’  H 


The  Laplace  transformation  formula  corresponding  to  the  flow 
equation  of  the  solenoid  valve  is: 

Q(s)=KeE(s)-KpHPH(s)  (20) 

3.  Equation  of  the  Thrust  Chamber  Head’s  Containing  Chamber 

When  the  pressure  of  the  thrust  chamber  head's  containing 

2 

chamber,  based  on  design  requirements  P„(t)  <  20kg/cn  ,  is  com- 

n 

pared  with  the  common  hydraulic  system,  the  containing  chamber's 
pressure  is  lower.  In  the  containing  chamber  equation,  it  is  not 
necessary  to  consider  the  containing  chamber  case  and  the  deform¬ 
ation  of  the  air  caused  by  the  pressure  in  the  containing  chamber. 
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Fig.  6#  The  hydraulic  containing  chamber 

Because  the  spring  has  prestress  F  , 

injector's  needle  valve  opens  and  moves, 

PH(t)=PH  •  y<t)=0-  Therefore, 
o 

ARdy(t)=Q(t)dt 

When  dy(t)  >  0,  Q(t)=Q1(t);  when  dy(t)  < 

y  (t)  =  —  Jo(t)dt 

The  corresponding  Laplace,  trafts format ion  formula  is: 

Y<«>-  SJpl  Q(s)  (22) 

4.  The  Equation  of  the  Injector; s  Needle  Valve  Movement 

Fig.  7  shows  the  stress  conditions  of  the  injector's  needle 
valve. 


Fig.  7  Stress  chart  of  the  injector's  valve. 


when  F  the 

that  is,  when 

0,  Q(t)=Q2(t). 

(21) 
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F 

O 

K 


y  (t) 

dy (t)/dt 


APH(ti 


pc<t) 


is  the  spring's  prestress  [kg] 

is  the  spring's  rigidity  [kg/mm] 

is  the  needle  valve's  displacement  [mm] 

is  a  certain  moment  's  movement  speed  of  the  needle 
valve  [mm/s] 

is  the  hydraulic  pressure  when  the  needle  valve  begins 
to  rise  [kg/mm2] 

2 

rs  the  increment  of  the  hydraulic  pressure  [kg/mm  ] 

2 

is  the  operating  area  of  the  hydraulic  pressure  [mm  ] 

is  the  damping  coefficient  of  the  needle  valve's 
movement  [kg* s/mm] 

is  the  operating  area  of  the  gas  towards  the  needle 
valve  [mm] 

2 

is  the  gas  pressure  [kg/mm  ] 

2 

is  the  needle  valve's  mass  [kg-s  /mm] 


The  equation  injector  needle  valve's  movement  is  (see  refer¬ 
ence  [  1  ] )  : 

2 

m  d  +B  ^+Kytt)*A  P  (t)+A.4Pu(t)  (23) 


The  corresponding  Laplace  transformation  formula  is: 


(mS  +BS+K)Y(s)=AcPc(s)+AHPH(s) 


5.  The  Equation  of  the  Injector's  Flow 


G  (t)=X  A 
o  o  o 


oxidizing  agent 


G^(t)=>t^Aj  2g  Yj  Pjs  incendiary  agent 


Many  tests  have  proven  that  when  the  displacement  of  the 
injector's  needle  valve  is  relatively  small,  the  changes  of  the 
flow  and  displacement  form  a  linear  relationship.  The  two  flow 
equations  can  be  expressed  as: 

GQ(t)«  ?0y(t) 

Gf(t)»  rfy(t) 
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(25) 

(26) 


In  the  formula,  coefficients  "f  and  can  be  determined  by 

'  o _  f 

flow  tests  of  the  injector. 

The  corresponding  Laplace  transformation  formula  of  the  flow 
equation  is: 

Go<S)=  V<S>  (27) 

Gf  (S)“  *fY(S)  (28) 

6.  The  Combustion  Chamber  Equation 

Basic  hypotheses: 

1)  Combination  delay  %  is  constant  and  it  is  considered  that 
passing  TT  time  after  the  propellant  is  sprayed  into  the  combus¬ 
tion  chamber,  the  gas  is  instantaneously  produced. 

2)  The  gas  is  an  ideal  gas. 

3)  When  comparing  the  flow  changing  speed  and  pressure  wave 
propagation  speed,  it  is  consiered  that  the  pressure  wave  propa¬ 
gation  speed  is  infinitely  great  and  the  pressure  of  the  combus¬ 
tion  chamber  at  each  location  is  uniform  at  any  given  moment. 

We  can  obtain  the  combustion  chamber  equation  (see  reference 
[1])  based  on  the  above  hypotheses. 

V  dp  (t).A. 

-af-  /P0<t>-G0Ct-t>«3£<t-T>  (29) 

In  the  formula, 

3 

v  is  the  volume  of  the  combustion  chamber  [mm  ] 
c 

R  is  the  constant  of  the  ideal  gas  [mm/°K] 

Tc  is  the  gas  temperature  [°K] 

2 

is  the  area  of  the  nozzle  throat  [mm  ] 

&  is  the  comprehensive  parameter  of  the  combustion 
chamber  [s] 
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The  corresponding  Laplace  transformation  formula  is: 


7. 


gf-S+^-  Pc(s)=e“rs[Go(s)+Gf  (s)  ]  (30) 

The  Equation  of  the  Feedback  Pressure  Sensor  and  Circuit 
Operational  Amplifier 


The  pressure  sensor  used  for  this  engine  system  is  a  resist¬ 
ance  type  pressure  sensor.  We  can  consider  that  it  is  an 
amplification  link  for  this  type  of  sensor.  Its  equation  is: 

Vf  (t)-KfPc(t)  (31) 

The  circuit  operational  amplif ier  equation  is : 

E(t)=Vc(t)-Vf  (t)  (32) 

The  corresponding  Laplace  transformation  formulas  are: 

Vf (s)=KfPc(s)  (33) 

E(s)=Vc(s)-Vf  (s)  (34) 

III.  The  Transfer  Function  of  the  Engine  Control  System 

We  synthesize  formulas  (20),  (22),  (24),  (27),  (28),  (30), 

(33)  and  (34)  and  use  a  block  diagram  to  indicate  the  relationship 
between  each  variab le. 


Fig.  8  Block  diagram  of  the  control  system. 
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In  the  diagram,  j=vc/RTc,  q=At46?  . 

Because  the  engine  uses  a  bipropellant  natural  propellant, 
ignition  delay  time  is  relatively  small.  Generally,  “T  K.  Sms 
which  can  be  considered  e  =1. 


Volume  v  of  small  engine  combustion  chambers  is  relatively 

small,  temperature  Tc  of  the  gas  is  relatively  high  and 

RT  v  .  It  can  be  considered  that  j=v  /RT  =0 . 

C  "  C  C  C 


The  block  diagram  is  simplified  into  the  following  form: 

K-K,V% 


W(s)  * 


e  c 
"5 - T 

ms  + (B+A^/Kp  )S+K-KcAc 
H 


(35) 


In  the  formula 


Kc»  i(Yo+Vf)=  (fo+ff) 


(36) 


K  is  called  the  combustion  chamber's  gain  [kg/mm  ]. 
c 


vt  <  s »  _£  <  s 


vf<  S) 


r 

W  iS> 

1 

L 

J 

> 

Kf 

P,  <  S  r 


Fig.  9  Block  diagram  of  simplification. 

The  closed-cycle  transfer  function  of  the  control  system  is: 


VcV*fH 


*  (S)  *ms5.  p  I^VycW1!. 

H  H 


(37) 


m 

B 


is  the  mass  of  the  injector's  needle  valve  [kg*svmm] 

is  the  damping  coefficient  of  the  injector's  needle  valve 
l  kg*  s/ram] 


A^  is  the  operating  area  of  the  needle  valve  [mm  ] 


V 
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is 

the 

solenoid  valve's 

pressure 

coefficient 

[  mm V  s • kg ] 

e 

is 

the 

solenoid  valve's 

flow  gain 

[mm3/s* V] 

is 

the 

spring's  rigidity 

[kg/mm] 

c 

is 

the 

combustion  chamber's  gain 

[kg/mm  ] 

i 

c 

is 

the 

operating  area  of 

the  needle  valve's 

gas  [mm3] 

■f 

is 

the 

gain  of  the  feedback  pressure  sensor 

[mm2V/kg] 

When  calculating  the  transfer  function,  the  B,  P„.  (when  t=0, 

ht 

P  =P  .)  and  K  can  be  determined  by  the  test  method  based  on  the 
Ht  HO  C 

concrete  structure  and  operating  parameters  of  the  engine. 

IV.  Analysis  of  the  Performance  of  the  Engine  Control  System 
1.  Analysis  of  the  Stability  of  the  Control  System 


The  characteristic  equation  of  the  control  system's  transfer 
function  is: 

mS2+  (B+A^/Kp  )  S+K-KcAc+KeKcKf  Ajj/Kp  *0 
H  H 

Letting  ~ 

*  -B+ Vkph’  V  =K-VV Vck£ VKPH 

The  root  of  the  characteristic  equation  is: 


S, 


—  A  +  ^  A1  — 4m» 
2m 


(38) 


—  A— s/A1  — 4  mv 
2m  ■ 


(39) 


Because  A  ^  0 ,  to  discuss  the  stability  of  the  control 
system  it  is  only  necessary  to  discuss  root  S^  of  the  character¬ 
istic  equation.  We  can  see  from  root  S^  of  the  characteristic 
equation  that: 

1)  When  v  <  0,  then  we  have  the  real  root  of  S^^O,  the 
output's  transition  component  has  monotonic  increases  based  on  the 
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exponential  law  and  the  system  is  unstable.  However,  even  if  the 
plan  of  this  engine  control  system  appears  unstable,  the  combus¬ 
tion  chamber's  pressure  will  not  be  able  to  have  monotonic 
increases  or  decreases.  Because  the  control  system  is  a  closed- 
cycle  system  with  feedback,  when  the  combustion  chamber's 
pressure  exceeds  a  certain  value,  this  can  cause  the  corresponding 
solenoid  valve  to  operate.  This  type  of  combustion  chamber  pres¬ 
sure  can  operate  according  to  certain  frequency  and  amplitude 
oscillation  mode.  The  exit  and  entrance  solenoid  valves  operate 
alternately. 

2 

2)  When  the  real  root  of  v>0,  7C  ‘-4mv^0  and  and  the 

control  system  is  stable. 

2 

3)  When  v^O,  X  -4mv^0,  and  S2  are  the  conjugate  complex 
roots  of  the  characteristic  equation  and  the  transition  component 
of  the  output  function  appears  as  an  oscillation  mode. 

In  order  to  further  display  that  these  parameters  are  major 
influencing  factors  of  the  dynamic  state  performance  of  the 
engine  control  system,  we  used  the  parameters  of  the  engine  to 
carry  out  concrete  calculations. 

Original  parameters: 


(1)  *  *  ( 

2)  *  *  |(3)  *  & 

!  (4)  *  « 

m  1  kf-i'/mm 

i.:xur* 

( 6 ) 

Am  i 

irxxo« 

(7)*W*Wld8®W 

Am  1  mm1 

2.SXJ0* 

(8 )  *(*!!£** 

B  i  kf-s/om 

M 

K  1  Icf/om 

m 
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4 


i 


(10)  *  fc  Clip  4  (12)*  «  |  (13)  *  * 


(14)  fcercxaiw 

A 

mm* 

i  0.64 

(15)  *e*;*ft*a 

i 

0.60 

( 16 ) 

P  Hi 

kg/mm1 

SXIO'1 

(17) 

P. 

|  kg/mtn’ 

i  18X10*’ 

(18)  *#!jr)*XfWR* 

\  f 

H* 

|  SO 

(19) 

T‘ 

i  s 

6X10-* 

(20) 

K, 

,  s/v 

1.2X10-’ 

(21)  *mMi8.EE 

E, 

V 

I  2 

(22)  57jf»S»*a 

K, 

i  V'3ia’/kg 

55.56 

Table  1 


Key:  (1)  Name;  (2)  Symbol;  (3)  Unit;  (4)  Numerical 
value;  (5)  Mass  of  the  injector's  needle  valve; 

(6)  Operating  area  of  the  needle  valve's  hydraulic 
pressure;  (7)  Operating  area  of  the  needle  valve's  gas; 
(8)  Damping  coefficient  of  the  needle  valve;  (9)  Rigid¬ 
ity  of  the  spring;  (10)  Name;  (11)  Symbol;  (12)  Unit; 
(13)  Numerical  value;  (14)  Flow  area  of  the  solenoid 
valve;  (15)  Flow  coefficient  of  the  solenoid  valve; 

(16)  Hydraulic  pressure  of  the  injector's  needle  valve; 

(17)  Entrance  pressure  of  the  solenoid  valve;  (18)  The 
circuit  pulse's  operating  frequency;  (19)  Pulse  width; 
(20)  Circuit  gain;  (21)  Voltage  of  the  deviation  signal; 
(22)  Gain  of  the  pressure  sensor. 


During  cold  tests,  we  used  water  as  the  hydraulic  pressure 
working  substance.  We  can  calculate  situations  with  rising 
thrust  based  on  the  data  given  in  Table  1: 

Ko -fKmAyj^Jl  -6.78x  lO'CmmYs-V-kg'''] 

-y/p'-PH,  -0.32[kg/mnr’]1/! 

K*“\  3-1  *=0.32CV-mrn/kg,/'] 

K'~KqK!  **2. 17  x  lO’Cmm’/s*  V] 

/C,w-/C0/C;-2.17x  10’[mm*/s.kg3 
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4 


If  the  injector  is  designed  well,  when  needle  valve  dis¬ 
placement  h  ^ 0.25mm,  the  change  of  the  injector's  flow  and 

ITlclX 

needle  valve's  displacement  is  close  to  a  linear  relationship. 
The  combustion  chamber's  gain  Kc  can  be  expressed  as: 


Pc 


Kc=  h 


max 


max 


-2  2  3 

If  Pc  ^  =10x10  kg/mm  and  h  =0 . 25mm,  then  K  =0.40  kg/mm  . 
max  max  c 

We  can  calculate  the  v  and  71  values  in  the  characteristic 
equation  based  on  the  above  mentioned  parameters  and  determine 
the  stability  of  the  engine  control  system.  Moreover,  we  can 
also  see  the  influence  of  each  parameter  on  the  transfer  function 
and  further  analyze  the  major  influencing  factors  on  the  dynamic 
state  performance  of  the  engine  control  system. 

r 

v  =  K-K.At+K.KcK,AH/KPli  =120- 100  +  355. 58  x  10*  =  35.578  x  lC’fcg 'm 
A-B+^/K- ,„=*3. 6  +  11. 797  x  10!«11.83x  lO’kg-s/'mn 
A* -4mv- 139.9489  x  10‘-1. 71-139. 487x  lO'Ckg-s/mm)* 

2 

Calculations  show  that  when  v  ^  0  and  71  -4mv  ^  0 ,  the  root 
of  the  transfer  function  characteristic  equation  is  a  negative 
real  root.  Therefore,  the  engine  control  system  is  stable. 

2.  Analysis  of  the  Transition  Time  of  the  Engine  Control 
System 


Because  the  mass  of  the  engine  injector's  needle  valve  is 

2 

very  small,  when  m  B+A^/Kp  ,  we  can  delete  the  second  degree 

**  *  tl 

term.  The  combustion  chamber  pressure  changes  measured  by  heat 
tests  have  monotonic  rises  and  decreases  based  on  the  exponent¬ 
ial  curve.  Therefore,  the  engine  control  system  can  be  simpli¬ 
fied  into  a  first  order  system  and  the  transfer  function  of 
equation  (37)  changes  into: 
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40) 


K.K€AhI K?„-v  _  g 
CA/OS+7  r5  +  1 


$=*K.K'Ah!K  ,h-v  (ftJ-il&.-k'&WL) 

T  =*i/v  (2) 


(1) 


Key:  (1)  Called  the  amplification  coefficient; 
(2)  Called  the  time  constant. 


Based  on  analysis  of  the  first  order  link,  when  combustion 

chamber  pressure  p  reaches  to  95%  of  the  specified  value,  its 

transition  time  t  =3T  and  therefore: 

P 


t .  3<B+ Vkp„  i 

p  K-*VVKeKcvH/KI 


(40) 


'H 


Because  the  range  of  change  of  each  parameter  is  not  large, 
we  can  see  from  the  previous  numerical  calculations  that: 

b«ah/kph'  l*-Ve'«WfV\ 


Relational  formula  (43)  of  transition  time  t  can  be  sim- 

P 

plified  into: 


3  A 


t- 


H 


P  K  K  K_ 
c  e  c  f 


(41) 


We  can  see  from  formula  (41)  that  the  transition  time  is 
mainly  determined  by  the  injector  needle  valve's  hydraulic  pres¬ 
sure  operating  area  A^  and  the  solenoid  valve's  flow  gain  Ke* 

In  order  to  reduce  transition  time  t  ,  we  reduce  as  much  as 

P 

possible  the  injector  needle  valve's  hydraulic  pressure  operating 
area  ,  appropriately  enlarge  the  solenoid  valve's  flow  area  A, 
increase  the  pulse  width  or  raise  the  pulse  operating  frequency. 
All  of  these  are  effective  measures  for  reducing  the  transition 
time. 

3.  Analysis  of  the  Steady  State  Errors  of  the  Engine  System 
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We  can  obtain  the  Laplace  transformation  expression  of 
deviation  signal  E(t)  from  the  block  diagram  in  Fig.  9  of  the 
engine  control  system: 


E  (s)  = 


1+W(s) K,Vc (s) 


In  the  formula 


W  (s )  = — 
ms 


VcV^H 

2+(B+A^/Kp  )S+K-KoA= 


Control  voltage  V  (t)  of  the  engine  control  system  introduced 
in  this  article  uses  the  jump  function  form  for  input.  There¬ 
fore,  we  will  only  discuss  the  positional  errors  of  the  jump 
input  signals  here  but  we  will  not  discuss  the  speed  errors  and 
acceleration  errors. 

V  (s) =Ky 
c  S 


In  the  formula,  is  the  amplification  coefficient  of  the  con¬ 
trol  voltage  [V] . 


We  can  obtain  the  general  expression  of  the  steady  state 
error  by  using  the  last  value  theorem  of  the  Laplace  transforma¬ 
tion. 

1 (t) =LimSE (s) =  1+W(0) K 
s  *^o  t 

By  substituting  in  W(0) ,  we  can  obtain  the  steady  state  error 
expression  of  the  engine  control  system: 

K-KA 

K-KrA.+K.K'K,AH/k, 7 


In  the  same  way, Jk-K  A  K  K  K.A  /K_  . 

c  c  e  c  1  H 

The  steady  state  error  expression  can  be  simplified  into  the 
following  form: 
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K-K  A 

1(t)=KeKcKfVKPHKv 


From  formulas  (9),  (10),  (11)  and  (12),  we  obtain: 


K  ”  E  (po"PHt) 
tt  r 


Then,  the  steady  state  error  expression  is: 


(K-K  A  ) E , 

l(t)=  — I — K 


WnVHtv 


(42) 


Enlargement  of  spring  rigidity  K  and  reduction  of  the  injector 

needle  valve's  hydraulic  pressure  operating  area  A  can  cause 

H 

steady  state  error  l(t)  to  increase. 

When  the  regulation  range  of  the  thrust  is  large  and  deviation 
signal  Efc  is  large,  then  the  steady  state  error  is  large.  The 
steady  state  errors  are  larger  when  there  is  high  thrust  (P  is 
large)  regulation  than  when  there  is  low  thrust  (P  is  small) 

nu 

regulation. 

The  situations  discussed  above  are  concerned  with  rising 

thrust.  As  for  situations  with  descending  thrust,  the  solenoid 

valve's  flow  gair  K  and  the  solenoid  valve's  pressure  coefficient 

Kp  are  substituted  into  relational  formulas  (15) ,  (16) ,  (17)  and 
H 

(18),  that  is,  by  using  the  exit  solenoid  valve's  parameters  and 
other  invariant  parameters,  we  can  obtain  the  same  analytical 
results. 

We  should  also  point  out  that  in  the  analysis  of  the  regulation 
performance  of  the  engine  control  system,  we  did  not  take  into 
consideration  the  steady  state  errors  of  the  components  of  the 
engine  control  system  which  are  caused  by  the  structural  parametric 
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changes  which  influence  their  sensitivity.  For  example,  when 
the  interval  of  the  solenoid  valve's  needle  valve  movements 
enlarge,  the  magnetic  resistance  enlarges  and  the  leading  edge 
of  the  open  needle  valve  enlarges,  these  cause  the  steady  state 
errors  to  increase.  In  the  same  way,  when  the  solenoid  valve's 
entrance  pressure  pQ  enlarges  and  the  starting  force  increases, 
this  can  also  cause  the  leading  edge  of  the  open  needle  valve  to 
enlarge  and  influence  the  dynamic  and  steady  state  performances. 
Furthermore,  the  nonlinearity  etc.  of  the  feedback  pressure 
sensor  can  also  influence  the  regulation  performance  of  the 
engine  control  system  but  these  factors  were  not  included  in  the 
above  analysis. 
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ELECTRONIC  SERVO  VALVE 
by  Cai  Yongnian 

Abstract 

This  paper  presents  a  new  type  of  electronic  servo  valve 
(ESV)  which  is  composed  of  a  pulse  modulator,  electro-mechanical 
actuator  and  variable  section  area  valves.  The  operating  prin¬ 
ciple  and  transfer  function  for  ESV  are  also  described.  At  the 
same  time,  we  also  give  test  data  on  ESV. 

I.  Introduction 

ESV  are  composed  of  servo  circuits  made  up  of  integrated 
circuits,  motors  and  variable  section  area  valves.  The  feedback 
effects  of  mounting  a  position  sensor  on  a  valve  rod  causes  this 
type  of  valve  to  have  relatively  high  static  precision  and 
relatively  fast  actuation  time.  With  the  use  of  millivolt  or 
volt  level  command  voltage,  it  can  control  the  flow  of  the  medium 
within  a  very  wide  range.  The  relationship  between  the  command 
voltage  and  flow  can  be  linear  and  it  can  also  be  nonlinear.  The 
servo  circuit  is  installed  in  a  case.  The  motor  and  variable 
section  area  valve  combine  into  an  organic  whole.  The  servo 
circuit  is  connected  by  the  cable  and  motor.  The  entire  structure 
of  the  ESV  is  shown  in  Fig.  1. 


Fig.  1  Entire  structure  of  the  ESV. 
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Because  this  type  of  valve  does  not  require  a  supplementary 
air  source  and  hydraulic  source  in  the  process  of  controlling 
the  flow,  this  type  of  valve  is  especially  applicable  in  situa¬ 
tions  which  do  not  have  an  air  source  and  hydraulic  source.  At 
the  same  time,  this  valve  is  also  applicable  for  control  of 
poisonous  mediums  and  situations  which  do  not  allow  people  to 
work.  Therefore,  this  valve  possesses  relatively  strong  struct¬ 
ural  adapted)  ility. 

One  of  the  specific  applications  of  the  ESV  is  its  use  to 
change  the  thrust  of  a  spacecraft's  liquid  rocket  engine  so  as  to 
cause  the  spacecraft  to  be  able  to  complete  soft  landing,  vari¬ 
able  orbit  flight  and  other  complex  space  navigation  tasks  on 
planets.  Figs.  (2) ,  (3)  and  (4)  show  the  situations  of  the 
actual  test  runs.  Fig.  2  is  the  situation  when  the  thrust  is 
relatively  large  and  because  of  this  the  corresponding  flame  is 
relatively  long.  Figs.  (3)  and  (4)  are  the  situations  when  the 
thrust  is  intermediate  and  relatively  small  and  because  of  this 
the  corresponding  flames  are  relatively  short.  In  Figs.  (2) , 

(3)  and  (4),  the  upper  part  of  the  combustion  chamber  has  two 
similar  ESV.  One  controls  the  oxident  of  the  engine  and  one  con¬ 
trols  the  fuel  of  the  engine. 


P-  C75<)97 


ar> 


P-67S<Wf 

— r~r ^  i  . 


L 


L. 


i 


P-67S9?? 


am 


Figs.  2,  3  and  4. 


II.  Fundamental  Principle 

The  ESV  is  composed  of  a  servo  circuit,  motor  and  variable 
section  area  valves.  The  position  sensor  is  installed  on  the 
end  of  the  valve  rod  and  uses  the  feedback  effects  of  the  posi 
tion  sensor  to  cause  this  type  of  valve  to  have  general  closed 
cycle  system  performance. 


Fig.  5  Chart  of  the  principle  of  the  ESV  system, 
(see  next  page  for  key) 
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Fig.  5 


Key:  (1)  Command  voltage;  (2)  Voltage  amplifier; 

(3)  Modulator;  (4)  Power  amplifier;  (5)  Torque  motor; 

(6)  Oscillator;  (7)  Displacement  sensor;  (8)  Compen¬ 
sator;  (9)  Demodulator;  (10)  Rolling  lead  screw; 

(11)  Variable  section  area  venturi. 

Dotted  line  frame  1  in  Fig.  5  is  the  servo  circuit.  It 
includes  the  voltage  comparator,  modulator,  power  amplifier, 
oscillator,  demodulator,  compensator  etc.  The  major  function  of 
the  servo  circuit  is  to  amplify  the  very  weak  error  signals  and 
transform  them  into  pulse  modulated  width  signals  which  can  drive 
the  motor.  The  frequency  of  the  pulse  is  determined  by  the  fre¬ 
quency  of  the  oscillator.  The  width  of  the  pulse  is  determined 
by  the  size  of  the  error  signal.  The  voltage  amplifier,  modula¬ 
tor  and  power  amplifier  are  positive  loops  of  the  closed-cycle 
system.  The  demodulator  and  compensator  are  feedback  loops  of 
the  closed-cycle  system.  Broken  line  frame  2  in  Fig.  5  is  the 
motor's  final  control  element.  It  is  composed  of  the  permanent 
magnet  type  torque  motor  and  the  rolling  lead  screw.  The  lead 
screw  transforms  the  rotating  movement  of  the  motor's  roll  into 
the  direct  movements  of  the  valve  rod.  The  guide  bar  of  the 
displacement  sensor  is  fastened  on  the  end  of  the  valve  rod.  The 
output  voltage  of  the  displacement  sensor  is  the  input  signals  of 
the  demodulator.  Broken  line  frame  3  of  Fig.  5  is  the  variable 
section  area  valve.  It  can  be  the  variable  section  area  venturi 
and  can  also  be  the  variable  section  area  current-limiting  ring. 
The  movable  valve  rod  extends  into  the  opening  of  the  venturi  or 
current-limiting  ring.  The  valve  rod  is  machined  into  a  certain 
taper  in  the  throat  area  near  the  venturi  and  the  smallest  aper¬ 
ture  of  the  current-limiting  ring.  Because  of  the  effect  of 
this  tapering,  when  the  valve  rod  moves  to  the  left  and  right, 
this  causes  the  flow  area  to  have  corresponding  changes  which 
changes  the  flow  quantity.  Structurally,  the  motor  and  valve 
sections  are  joined  by  the  mechanical  method.  Figs.  (6)  and  (7) 
are  detailed  blueprints  of  this  type  of  structure. 
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Fig.  5  Full  diagram  of  the  motor  mechanism  and  variable  section 
area  venturi. 
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Fig.  7  Full  diagram  of  the  motor  and  the  variable  section  area 
current-limiting  ring. 

The  structure  in  Fig.  6  uses  the  venturi  to  control  the  flow 
quantity.  The  structure  in  Fig.  7  uses  the  current-limiting 
ring  to  control  the  flow  quantity.  The  structures  of  the  two 
are  basically  the  same,  only  the  structures  of  the  valve  sections 
are  different.  The  external  dimensions  of  Fig.  7  are  smaller 
than  those  of  Fig.  6. 
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The  operating  principle  of  the  ESV  is  as  follows.  When  the 
system  is  in  a  holding  state,  the  error  signal  in  front  of  the 
voltage  amplifier  is  zero  (in  reality,  it  is  not  zero  but  very 
small) .  Therefore,  after  passing  the  modulator  and  power  ampli¬ 
fier,  the  pulse  signals  added  to  the  torque  motor  are  pulse 
signals  with  equal  widths.  Because  the  mean  value  of  the  pulse 
signals  with  equal  widths  is  zero,  the  torque  motor  does  not 
rotate.  Therefore,  the  valve  rod  is  also  in  the  static  state. 

At  this  time,  the  output  signals  of  the  displacement  sensor  are 
zero.  When  the  command  mechanism  gives  a  certain  command  volt¬ 
age  because  the  output  signals  of  the  displacment  sensor  are 
zero,  error  signals  are  produced.  The  error  signals  are  ampli¬ 
fied  and  modulated  by  the  positive  loop  which  causes  the  error 
signals  to  change  into  corresponding  pulse  signals  with  unequal 
widths.  Because  the  mean  voltage  of  the  pulse  signals  with 
unequal  widths  is  zero,  the  motor  rotates.  The  rolling  lead 
screw  rotates  and  changes  the  movements  of  the  motor  into  direct 
movement  of  the  valve  rod.  Following  the  movement  of  the  valve 
rod,  the  displacement  sensor  then  gives  the  corresponding  electric 
signals.  With  the  effects  of  the  feedback  loop,  these  electric 
signals  are  transformed  into  corresponding  position  feedback 
signals  and  they  are  compared  to  the  command  voltage.  After 
comparison,  the  obtained  error  signals  are  again  regulated  ac¬ 
cording  to  the  above  mentioned  process.  When  the  regulation 
process  is  completed,  the  error  signals  are  zero,  the  valve  rod 
stops  moving  and  each  position  of  the  valve  rod  corresponds  to 
each  command  voltage.  The  relationship  between  them  is  zero. 

III.  The  Transfer  Function  of  the  ESV 

Because  the  response  time  of  an  electronic  package  is  much 
faster  than  that  of  the  motor  and  valve,  when  the  time  constant 
of  the  electronic  package  is  omitted,  we  simplify  Fig.  5  into 
Fig.  8.  We  deduce  its  transfer  function  based  on  Fig.  8 
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Pig.  3,  Block  diagram  of  the  ESV. 


We  can  see  from  Pig.  8  that  the  terminal  voltage  of  the  two 
terminals  of  the  motor's  armature  can  be  expressed  by  the  fol¬ 
lowing  formula: 

r»)KAK, 


The  differential  equation  of  the  armature  circuit  is  as  follows: 


E„=  L 


'*■  dt 


+R*i4+eb 


Torque  T  produced  when  the  armature  was  under  a  corresponding 
speed  is  in  direct  ratio  to  the  armature's  electric  current.  Its 
relational  formula  is  as  follows: 

T=KTi* 

Counter  electromotive  force  in  the  armature  loop  can  be  ex¬ 
pressed  by  the  following  formula: 

eb=Kb  if 

The  armature  rotates  at  a  corresponding  speed.  At  the  same  time, 
the  produced  torque  will  satisfy  Newton's  second  law.  The 
obtained  equation  is  as  follows: 

T“J  ^  +£  +Vtl 

The  transfer  effect  of  the  lead  screw  causes  the  friction 
coefficient  of  the  valve  section  converted  on  the  motor's  axis 
to  change  very  little.  Thus,  the  f  in  the  above  formula  is 
actually  the  viscous  friction  coefficient  of  the  motor's  axis. 
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Tc  is  the  friction  torque  converted  on  the  motor's  axis  caused 
by  the  dynamic  seal  of  the  two  terminals  of  the  valve  rod. 


(When)  Te-+Mc 

(When)  Tc=-Mc 


in  the  above  formula  is  the  friction  torque  caused  by  the 
coulomb  friction.  is  the  load  torque  of  the  fluid  hydraulic 

power  converted  on  the  motor's  axis.  The  transfer  effect  of  the 
rolling  lead  screw  causes  the  driving  torque  in  the  aforementioned 
formula  to  be  Tc+Tl>  In  order  to  simplify  the  calculation  and 
analysis  of  the  simplified  system,  we  omitted  the  effects  of 
Tc+Tl.  In  this  way,  we  caused  the  nonlinear  differential  equation 
to  be  simplified  into  a  linear  differential  equation.  Its  formula 
is  as  follows: 


T=J 


,2b  da 


+f 


dt 


dt~ 


Drawing  support  from  the  drive  effects  of  the  lead  screw,  the 
rotating  movements  of  the  motor's  armature  change  into  linear 
movements  of  the  valve  rod.  Its  relational  formula  is: 

jV©  -8 


After  Laplace  transformation  and  arrangement,  the  block  diagram 
after  the  Laplace  transformation  is  as  shown  in  Fig.  9  and  the 
transfer  function  corresponding  to  it  is  as  follows: 


H (S)  _ 


K. K_K_  -b 
A  P  T  27T 


S [LaJS*+ <Laf +RaJ) S+Raf +KTKb] +KAKpKTKFB 
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Fig.  9  Block  diagram  of  the  Laplace  transformation  of  the  ESV. 


The  root  locus  of  the  above  as  yet  uncompensated  transfer  function 
is  shown  in  Fig.  10.  After  KpB  and  are  selected,  when  K^, 

Kp  and  KT  are  relatively  low,  the  system  indicated  by  the  root 
locus  chart  is  stable.  However,  when  K^,  Kp  and  Kt  become  very 
large,  the  system  is  unstable  which  was  proven  by  tests. 


Fig.  10  Root  locus  chart  of  the  ESV. 

IV.  Selection  of  the  Unit 

The  comparison  amplifier  is  the  convergent  area  of  the  command 
signals,  position  feedback  signals  and  compensating  signals.  Its 
reliability  and  stability  will  directly  influence  the  performance 
of  the  system.  The  polarity  of  the  command  signals  and  the  pol¬ 
arity  of  the  feedback  signals  are  opposite  when  the  frequency  is 
relatively  low.  The  compensating  signals  are  used  to  improve  the 
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performance  of  the  system.  The  above  mentioned  comparator 
amplifier  uses  a  common  solid  state  unit  -  the  operational 
amplifier.  See  Fig.  11  for  its  principle.  Its  output  voltage 
Vq  can  be  indicated  by  the  following  expression: 


v  COm 


Fig.  11  Chart  of  the  principle  of  the  voltage  comparator. 

Key:  (1)  Command;  (2)  Feedback;  (3)  Compensator; 
(4)  Output. 


In  the  figure,  R^,  R2  and  R^  are  isolation  resistances. 
R^=R2=R3»  then  the  above  formula  changes  to: 


v= rv  +v  ,  f 

O  FB  CWr^ 


When 


The  comparison  amplifier  amplifies  the  error  signals  and  after¬ 
wards  sends  them  to  the  modulator.  The  modulator  transforms  the 
error  signals  into  equally  high  pulse  signals.  The  modulator 
also  uses  the  common  operational  amplifier.  The  input  terminal 
of  the  modulator  is  the  3um  total  of  the  error  signals  and 
triangular  waves.  Use  of  the  high  gain  and  limiting  effects  of 
the  operational  amplifier  causes  the  input  signals  of  the  modu¬ 
lator  to  transform  into  equally  high  rectangular  pulses  with 
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variable  widths.  The  chart  of  its  principle  is  shown  in  Figs. 
(12)  and  (13). 


Fig.  12  Chart  of  the  principle  of  the  modulator. 


Fig.  13  Chart  of  the  waveforms  of  the  modulator. 


If  e  is  the  output  voltage  of  the  comparison  amplifier,  efc 
is  the  triangular  wave’s  voltage,  the  triangular  wave  peak-peak 
voltage  amplitude  is  +  Ep  and  the  void  ratio  occupied  by  the 
modulator's  output  pulses  in  R,  then  R  can  be  expressed  as  follows 


R+0.5- 


In  the  formula,  T  is  the  pulse  period  and  t1  is  the  width  of 


the  modulated  pulse.  In  order  to  obtain  the  linear  relationship 
between  R  and  e  ,  the  linearity  of  the  triangular  waves  must  be 
very  good. 

The  power  amplifier  of  the  driving  motor  uses  a  switching 
amplifier.  The  output  power  of  this  type  of  circuit  can  be  very 
large  yet  the  generated  heat  of  the  output  power  tube  is  very 
small.  Its  principle  is  shown  in  Fig.  14. 


Fig.  14  Chart  of  the  principle  of  the  switching  amplifier. 


Its  operating  principle  is:  when  the  T^  and  T^  are  conducting, 
the  armature's  current  is  the  i  ,  shown  in  Fig.  14.  When  the  T.. 
and  T^  powe,r  tubes  are  conducting,  the  current  is  i^*  there¬ 
fore,  when  the  conducting  time  of  the  power  tubes  is  changed, 
this  can  cause  the  motor  to  have  positive  and  negative  rotation. 
In  one  pulse  period,  the  voltage  of  the  two  ends  of  the  armature 
is  E  in  the  t1  time  interval.  in  the  T-t.  time  interval,  the 
voltage  of  the  two  ends  of  the  armature  is  -E  .  Because  of  this, 
within  one  pulse  period,  the  differential  equation  of  the  arma¬ 
ture  circuit  can  be  expressed  as  follows: 


di 

La  cTF-  +iaRa+eb=Ea  0  ^  ^  t' 

di 

L  +i  R  +e,  =  -E 

a  dt  a  a  b  a 


t  <  T 
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It  is  well  known  that  the  current  solution  of  the  above  equation 
is  exponentially  shaped.  Therefore,  comparison  is  complex.  In 
order  to  simplify  the  solution  of  the  equation,  let  us  first 
consider  the  item.  The  steady  state  formula  of  e^  is: 

e,  =K.  2  'ff  n 
b  b 

In  the  formula,  K.  is  the  reaction  coefficient  of  the  armature 

b 

which  is  determined  by  the  structure  of  the  motor  and  n  is  the 
rotating  speed  of  the  armature.  In  order  to  bring  the  advantages 
of  the  torque  motor  into  full  play,  its  operating  state  should 
approach  the  blocked  rotating  point  and  also  rotate  at  very  low 
rotating  speeds.  In  order  to  satisfy  the  100  millisecond  time 
interval  required  for  the  entife  process  of  valve  movement,  the 
minimum  rotating  speed  of  the  motor  must  be  greater  than  12 
rotations/second.  However,  the  integration  time  of  the  afore¬ 
mentioned  equation's  solution  of  the  current  (i.e.  one  pulse 
period)  is  very  small.  The  pulse  period  T  selected  for  this 
valve  is  less  than  0.5  milliseconds .  This  time  interval  is 
much  smaller  than  the  time  required  for  one  turn  of  the  armature 
(•j-j  second  =  83  milliseconds) .  Therefore,  it  is  reasonable  to 
consider  that  within  the  integration  time,  the  rotating  speed  of 
the  motor  is  constant  and  invariant.  Thus,  it  is  assumed  that 
e^  does  not  change  within  one  pulse  period.  Secondly,  in  design, 

we  soucht  to  have  the  armature's  resistance  R  be  as  low  as 

If  a 

possible.  If  we  used  mean  current  I  =  —  ix_(t)dt  to  represent 

a  t  j  a 

i  ,  the  product  of  R  I  is  very  small  and  thus  the  proportion 

cl  a  a 

occupied  by  it  in  V  is  also  very  small.  If  we  use  mean  voltage 

V  to  represent  the  sum  of  the  above  two  items,  then  we  can 

approximately  consider  that  V  =i  R  +e.  =1  R  +e.  is  constant  and 

a  a  a  b  a  a  b 

invariant.  The  approximate  solution  of  the  above  obtained  cor¬ 
responding  armature  current  is  shown  in  the  following  formula: 

E  -V 

ia  (t)  =ia  (°)  +  — — -t  0$  t  ^  tx 

a 

E  +V_ 

ia(t)=ia(tl)~  ~I — “  “(t-ti)  t1  t  T 

cl 
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The  current  changes  appearing  on  the  two  ends  of  the  armature 
are  indicated  by  the  following  formula: 


E  T 
a 


(1  -P2) 


In  the  formula,  /O  =VC  .  When  V  =0,  then  /O*0.  At  this  time 

V-  L 

Cmax 

the  current  changes  of  the  two  ends  of  the  armature  are  maximum 
values : 

E  T 
_  a 

amax  2L 

c 

We  can  see  from  the  above  that  when  the  motor  does  not  rotate, 
there  is  an  alternating  current  flowing  in  the  armature.  The 
generation  of  heat  is  not  advantageous.  Yet,  because  there  are 
continuous  positive  and  negative  direction  currents  in  the  arma¬ 
ture,  the  armature  is  as  if  in  a  "dynamic  lubrication"  state. 
This  type  of  state  is  advantageous  to  restarting  the  motor. 


The  displacement  sensor  is  the  main  sensitive  element  of  the 
feedback  loop.  In  order  to  increase  the  anti-interference  per¬ 
formance  of  the  system,  the  output  signals  of  the  displacement 
sensor  are  detected  by  the  phase  demodulating  method.  Naturally 
increasing  the  resolution  of  the  displacement  sensor's  output 
signals  and  decreasing  the  phase  shift  of  the  output  signals 
will  further  improve  the  performance  of  the  valve  system. 


In  order  to  reduce  the  action  time  and  lighten  the  weight  of 
the  mechanical  drive,  this  plan  uses  the  permanent  magnet 
torque  motor  to  drive  the  motor.  Because  this  type  of  motor  pos 
sesses  low  speed  and  high  torque  performance,  this  motor  avoids 
the  use  of  a  decelerator  and  thus  reduces  the  space.  Use  of  a 
combination  of  this  type  of  torque  motor  and  rolling  lead  screw 
can  cause  the  relatively  small  rotating  moment  to  transform  into 
relatively  large  thrust  or  pulling  force.  See  the  figures  in 
Pig.  6  and  the  right  half  of  Fig.  7  for  details  of  the  structure 
Number  1  in  Fig.  6  is  the  torque  motor,  and  numbers  7  and  15 
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are  rolling  lead  screws.  Number  36  in  Fig.  7  is  the  torque 
motor  and  symbol  37  is  the  rolling  lead  screw.  The  relationship 
between  thrust  F  on  the  lead  screw  and  torque  T  of  the  torque 
motor  is  expressed  by  the  following  formula: 


) 


In  the  formula,  K  is  the  empirical  coefficient  which  commonly 
uses  0.1  to  0.2  and  L  is  the  spin  moment,  rj  is  the  drive  power 
and  in  most  situations,  (K+i)  1.  Therefore,  F^  T  and 

relatively  small  torque  can  produce  relatively  large  thrust. 


In  this  system,  the  variable  section  area  valve  has  two  struct¬ 
ures.  One  is  the  variable  section  area  venturi.  Its  structure  can 
be  seen  in  Fig.  6.  The  other  is  the  variable  section  area  current- 
limiting  ring.  Its  structure  can  be  seen  in  Fig.  7.  In  a  certain 
position,  the  valve  rods  of  the  two  are  machined  into  a  certain 
taper.  In  the  movement  process  of  the  valve  rod,  due  to  the 
effects  of  this  taper,  the  circular  flow  area  continually  changes. 
The  effects  begin  from  the  surface  and  control  the  flow.  After 
deductions  and  arrangements,  the  relational  formula  of  distance  H 
of  the  valve  rc* ' s  movement  and  controlled  flow  quantity  Q  is: 

Q  =  li-n-H-  *ina(D  —  ti  sinacosa)  ^  2gY(P,—  £  .,) 

The  above  relational  formula  is  suitable  for  the  variable  section 

area  venturi  as  well  as  for  the  variable  section  area  current- 

limiting  ring.  When  it  is  a  variable  section  area  current- 

limiting  ring,  the  P  in  the  above  formula  should  change  to 

o  s 

P_.  P  is  the  fluid's  saturated  steo,n  pressure.  P0  is  the  down- 
stream  pressure  of  the  valve.  When  half  tapered  angle  a  of  the 
valve  rod  is  relatively  small,  the  Q-H  properties  of  the  variable 
section  area  venturi  are  linear. 

V.  Tests  and  Data 

This  valve  underwent  water  flow  tests,  medium  tests  and 
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ignition  test  runs.  The  systems  chart  of  the  water  flow  tests 
is  shown  in  Fig.  15. 


( 5 ) 


Fig.  15  Systems  chart  of  the  water  flow  tests. 

Key:  (1)  Air  bottle;  (2)  Water  tank;  (3)  Manual 
operated  valve;  (4)  Servo  circuit;  (5)  Electric 
source;  (6)  Command  voltage;  (7)  Flow  meter; 

(8)  Manual  operated  valve. 


The  entrance  pressure  is  indicated  by  P^.  The  exit  pressure  is 
indicated  by  There  are  three  flow  meters  in  the  system's 

downstream  and  their  aim  is  to  measure  the  flow  quantity  in  a 
very  wide  range.  Fig.  16  shows  the  degree  of  opening  character¬ 
istics  of  the  variable  section  area  venturi. 


Fig.  16  Opening  characteristics  of  the  variable  section  area 
venturi. 

Key:  (1)  Grams/second,  flow;  (2)  Casing:no.  5;  (3)  Rod 
number:  no.  5;  (4)  Opening  characteristics  of  the  var¬ 
iable  section  area  venturi;  (5)  Millimeters. 
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Fig.  17  is  the  steam  corrosion  characteristics  of  the  variable 
section  area  venturi.  The  vertical  coordinate  is  the  flow 
quantity  and  the  horizontal  coordinate  is  the  reverse  voltage 
of  the  downstream.  The  solid  lines  and  points  of  inflection  in 
the  figure  are  the  test  values.  The  broken  lines  are  elongated 
lines.  Each  line  in  the  figure  corresponds  to  the  degree  of 
opening  of  the  valve.  We  can  see  from  the  curve  that  the  gas 
corrosion  characteristics  of  this  type  of  venturi  increase  with 
the  decreases  of  the  degree  of  opening  of  the  valve. 


Fig.  17  Gas  corrosion  characteristics  of  the  variable  section 
area  venturi. 

Key:  (1)  Grams/second,  flow;  (2)  Casing:  no.  5;  (3)  Rod 

number:  no.  5;  (4)  P=30  kg/cm2;  (5)  Gas  corrosion 
characteristics  of  the  variable  section  area  venturi; 

(6)  Millimeters;  (7)  Kg/cm2. 

Fig.  18  shows  the  characteristics  of  the  degree  of  opening  of  the 
variable  section  area  current-limiting  ring.  The  vertical  coor¬ 
dinate  is  the  flow  quantity  and  the  horizontal  coordinate  is  the 
command  voltage.  When  comparing  the  test  processes,  it  should 
be  pointed  out  that  the  control  precision  of  the  variable  section 
area  venturi's  flow  quantity  is  relatively  high  and  the  pressure 
decrease  of  the  valve  is  relatively  large;  the  influence  of 
reverse  pressure  on  the  flow  quantity  of  the  variable  section 
area  current-limiting  ring  is  relatively  large  and  thus 
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fluctuation  is  very  large.  However,  the  pressure  decrease  of 
the  valve  can  be  very  small. 


Fig.  18  Characteristics  of  the  degree  of  opening  of  the  variable 
section  area  current-limiting  ring. 

Key:  (1)  Flow  quantity;  (2)  Grams  second;  (3)  Grams/ 
second;  (4)  Grams/second;  (5)  Servo  amplifier  number:2; 
(6)  Adjustable  valve  number:  5;  (7)  Characteristics  of 
the  degree  of  opening  of  the  variable  section  area 
current-limiting  ring;  (8)  Command  voltage. 


VI.  Conclusion 

After  the  first  stage  of  design  and  production,  tests  showed 
that  the  ESV  could  be  used  to  continuously  change  the  flow 
quantity  of  the  fluid.  By  using  it  we  can  compose  a  constant 
value  automatic  regulating  system  of  the  flow  quantity  and  a 
follow-up  system  of  the  flow  quantity.  After  the  ESV  is  added 
on  to  the  combustion  chamber,  we  can  compose  a  bipropeller  or 
single  propeller  variable  thrust  engine.  Its  test  conditions  can 
be  seen  in  Figs.  2,  3  and  4.  This  valve  is  also  suitable  for  use 
in  other  situations  which  require  continuous  control  of  the  flow 
quantity. 

Explanation  of  Symbols  Used  in  This  Article 
V  voltage 

K  amplification  coefficient 
R  resistance 
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T  torque,  period 

t  time  variable 

L  inductance 

J  mechanical  inertia 

S  Laplace  transformation  operator 
o  half  taper  angle 

u  flow  quantity  coefficient 

Q  flow  quantity 

P  pressure 

Y  specific  gravity 

g  gravitational  acceleration 
F  Force 


Explanation  of  Lower  Symbols 

C  command 
FB  feedback 
A  amplification 
P  power  amplification 
com  compensation 
a  armature  loop 
L  load 
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THE  SELECTION  OF  LAUNCH  TIME  FOR  NEAR-EARTH  SATELLITES  (SPACE¬ 
CRAFT)  WITH  MISSION  OF  VISIBLE  PHOTOGRAPH 

by  Fan  Jianfeng 
Abstract 

This  paper  presents  a  set  of  formulas  which  can  determine 
the  launch  time  of  the  visible  photographic  satellite.  If  these 
formulas  are  applied  to  determine  the  launch  time,  the  satellite 
will  have  the  maximum  photographic  latitude  or  a  given  area  will 
be  photographed  in  a  specific  time. 

I.  Preface 

When  satellites  equipped  with  visible  light  cameras  carry 
out  observation  photography  of  the  earth,  if  the  area  does  not 
have  good  light  or  it  is  a  dark  night,  then  the  satellite  will  be 
unable  to  take  good  photographs  or  the  film  will  come  out  blank. 
This  article  studies  the  laws  of  the  spot  below  the  satellite 
being  illuminated  by  solar  light  when  an  observation  satellite 
is  taking  visible  photographs  of  the  earth.  Based  on  these  laws, 
selection  of  the  appropriate  launch  time  can  satisfy  the  various 
requirements  of  visible  photography.  For  example,  obtaining  the 
maximum  photographic  area  or  photographing  specified  areas  at 
specific  local  times. 

II.  The  Coordinate  System 

In  selecting  the  eguatorial  coordinate  system,  the  position 
of  point  B  is  indicated  by  right  ascension  and  declination 
SP  .  The  right  ascension  is  measured  from  Spring  equinox  point 
X  along  the  equator  towards  the  east,  from  0°  to  360°.  The 
declination  is  measured  along  the  equatorial  circle,  from  0°  to 
+  90°.  From  the  equator  north  is  positive  and  towards  the  south 
is  negative  as  shown  in  Fig.  1. 
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Fig.  1 

Key:  (1)  North  pole;  (2)  Equator. 

III.  The  Photographic  Arc 

If  the  right  ascension  of  the  sun  S  is  a$  and  the  declin¬ 
ation  is  (£@>  ,  then  the  projection  of  the  demarcating  line  (not 

including  twilight  shadows)  of  the  brightness  and  darkness  of 
the  earth  on  a  celestial  sphere  is  great  circle  I  as  shown  in 
Fig.  2. 


Fig.  2 

We  call  this  great  circle  I  the  shadow  area  great  circle.  Its 
equation  is: 
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-ctg  P  cos  (a© -7L)  =tg<£© 


(1) 


If  the  included  angle  of  the  satellite's  orbital  surface 
and  the  surface  of  the  equator  is  i  and  the  angular  distance 
from  Spring  equinox  point  r  to  ascending  node  A  is  £2  ,  the 
projection  of  the  satellite's  orbit  on  the  celestial  body  is 
great  circle  II  as  shown  in  Fig.  2.  We  call  great  circle  II  the 
orbital  great  circle.  Its  equation  is: 

sin( Pi  -  £2  )ctg*P  =ctgi  (2) 


In  order  to  determine  photographic  arc  BUB'  (see  Fig.  2)  , 
we  must  find  the  positions  of  nodes  B  and  B'  of  great  circles  I 
and  II.  For  this  reason,  we  introduce  variable  ^  and  let: 

tg<£  =tg<P/cosA» 

Therefore,  formulas  (1)  and  (2)  can  be  written  as: 


.....  + _ tgjr _ 

ctg A©  cos  A.©  ctg<5<5> 


=  -1 


tg7l 


,  ,  =i 

siny  tgi 


By  solving  formulas  (1A)  and  (2A) ,  we  can  obtain: 

tqA  =  *3  4s  tgising  S2*A$L.  (3, 

y  tg<£©  tgicos^  +  sin*£© 

tg«P  =  i -  t|ico.(y-g) - ___  (31 

[1+tg  O©  tg  i+2tgQ®  tgisin(<a^  -&£]  ' 

Therefore,  we  obtain  the  coordinates  of  nodes  B  and  B'  of  great 
circles  I  and  II: 


V  tg  <3  ®  tg » cos  +  sin  a  * 


« , »  lg - ( _ tg>Cos(g®-fl) _ \ 

\  [1  +  tg:<3®  tg,»  +  2tg<3®  tgi  sin  (a®  — 


V-180*+2# 
9»,mm  ~<Pa 


Solar  coordinates  and  in  the  above  formulas  can  be 
looked  up  in  an  astronomical  calendar.  After  satellite  injectors, 
orbit  inclination  can  be  considered  invariant  and  ascending  node 
angular  distance  Q  is  composed  of  two  parts:  initial  value  ^o 


and  the  oblateness  of  the  earth  perturbation 


AQ. 


IV.  The  Relationship  of  Qo  and  the  Launch  Time 

The  injection  time  of  the  satellite  influences  Q  and  the 
ascent  orbit  of  the  satellite  can  usually  be  considered  to  be 
given.  Therefore,  as  soon  as  the  launch  time  is  given,  the  in¬ 
jection  time  is  also  determined  and  thus  the  launch  time  in¬ 
fluences  Q  . 

o 

If  Z  is  the  orbit  insertion  point,  V* is  the  latitude  of  the 
orbit  insertion  point  and  y is  the  orbit  insertion  direction 
(measured  from  north  to  east)  as  shown  in  Fig.  3. 


Fig,  3 

From  Fig.  3  we  can  obtain: 

Q  +  -(12h-t  )  (5A) 

O  2  Z 

In  the  formula,  12^  represents  12  hours,  t  represents  the  in- 

z 

jection  time  (local  time  of  the  orbit  insertion  point).  In  the 

formula,  has  the  following  relational  formula: 

z 

tg  A7L  =tg  V'  sinS9 

Z  Z 
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(5B) 


Therefore,  we  can  obtain: 

QQ=a,0-12h-tg  1(tgV'sin<J°  z)+tz 

V.  Local  Time  of  Each  Point  on  the  Photographic  Arc 

If  M  is  an  arbitrary  point  on  the  photographic  arc  and  its 
latitude  is  then  we  have  a  formula  analogous  to  formula 

(5A)  : 

£+  2MM=a,0-U2h-tM) 

In  the  formula,  t,.  is  the  local  time  of  point  M  and  A*'..  has  the 

M  M 

following  relational  formula: 

sin^Tl  M=ctgitg9PM 

Thus,  we  can  obtain: 

Q  =0,0  -12h-sin'1  (ctgitg  ¥>  M>  +tJ( 

By  using  formula  (5B) ,  we  finally  obtain: 

tM=sin  1(ctgitgSDM) -tg  1  ( tg<P  sin^  g)  +  A  £2+tz  (6) 

Because  orbit  inclination  i  only  depends  on  the  orbit  inser¬ 
tion  point's  latitude  and  the  orbit  insertion  direction,  we  can 
see  from  formula  (6)  that  the  local  time  of  each  point  of  the 
photographic  arc  is  related  to  the  injection  time,  orbit  inser¬ 
tion  direction  and  orbit  insertion  point's  latitude.  In  most 
situations,  the  orbit  insertion  point's  latitude  and  the  orbit 
insertion  latitude  are  often  limited  by  the  demands  of  each  area. 
Thus,  if  we  want  to  photograph  assigned  areas  in  specific  local 

times,  we  can  select  an  appropriate  injection  time  t  so  that 

z 

t„  in  formula  (6)  has  the  desired  value. 

M 

VI.  The  Maximum  Photographic  Area 

If  satellites  with  missions  of  visible  photograph  are  not 
used  to  observe  specific  areas  on  the  ground  but  carry  out  com¬ 
mon  observations  of  various  areas  on  the  ground,  then  it  is  hoped 


that  the  satellite  can  photograph  as  many  ground  areas  as 
possible. 

1.  The  Situation  of  Observation  Satellites  With  Missions 
of  Photograph  in  Short  Flights 


Because  the  photographic  film  carried  by  most  satellites  with 
missions  of  photograph  is  limited,  flight  time  does  not  exceed 
several  days.  Therefore,  when  analyzing  the  maximum  photographic 
latitude,  we  can  take  &Q  and  as  constants.  In  this  way,  we 

find  the  derivative  of  the  tgSPg  for  S2  from  formula  (3B)  and 
letting  the  derivative  be  zero,  we  can  obtain: 

sin2  (OU&-Q)  +  (tgitg6©  +ctgictgcSfi>  )  sin  -  Q)  +1=0 


By  solving  the  above  formula,  we  obtain: 

sin  («,$-<£)=  -tgitg  S(S>  (i^f  90°-/^/  )  (7A) 

or 

sin -ctgictg<f<9  (i^90°  !S@J)  (7B) 

Therefore,  there  are  not  many  flight  days,  we  can  eliminate  AQ. . 
use  to  replace  &  and  obtain: 


sin_,(  —  tgitg  3#) 

fl,,  =  a«+  dn-l(  —  tgitgi*)  — 180' 


(iOO*  —  |3®  | )  (8 A} 


or 


— ctgi'ctg<J*  )  ■  i 

O,t-.a*  +  *m-'(-ctgictg<5.)r180*  )  («>90* -  |<5.  | )  (8 B) 

Now  we  will  explain'  the  two  situations  corresponding  to  the 


two  solutions. 


By  substituting  the  solution  given  in  formula  (7A)  into 
formula  (3B) ,  we  can  obtain  maximum  photographic  latitude 


tg<P 


=  +tgi 
max  —  ’ 


At  this  time,  the  corresponding  photographic  latitude  range 
is:  [ i , — i ] .  This  type  of  situation  corresponds  ^o  the  situation 
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wherein  the  inclination  of  the  orbit  surface  is  smaller  than 
the  inclination  of  the  shadowed  great  circle  surface.  Ac  tais 
time,  the  nodes  of  the  two  great  circles  fall  on  the  apex  areas 
of  the  orbit's  great  circle  (i.e.  the  highest  and  lowest  latitude 
areas  of  the  orbit's  great  circle). 

By  substituting  the  solution  given  in  formula  (7B)  into 
formula  (3B)  ,  we  can  obtain  maximum  photographic  latitude  ¥  : 

±ct<3$® 

The  corresponding  photographic  latitude  range  is: 

<p  [i,  -90°  +  &e  ]  &e>  >  o 

<P  [90°  +  8<s>  ,  -i]  <  o 

This  type  of  situation  corresponds  to  the  situation  wherein 
the  inclination  of  the  satellite's  orbit  surface  is  larger  than 
the  inclination  of  the  shadowed  great  circle  surface.  At  this 
time,  the  nodes  of  the  two  great  circles  fall  on  the  apex  areas 
of  the  shadowed  great  circle  (i.e.  the  highest  and  lowest  latit¬ 
ude  areas  of  the  orbit's  great  circle). 

The  following  deductions  will  obtain  the  maximum  photographic 
area  and  the  proper  injection  time  of  the  satellite. 

We  can  see  from  formula  (5B)  that  when  the  launch  date,  orbit 

insertion  direction  and  orbit  insertion  point  of  the  satellite 

are  given,  SBn  is  determined  by  injection  time  t  .  Therefore, 
o  z 

by  properly  selecting  the  injection  time  (or  corresponding  launch 

time)  we  can  cause  5^  to  possess  our  desired  value.  For  ex~ 

o 

ample,  by  causing  to  satisfy  formula  (8A)  or  formula  (8B) , 

we  can  obtain  the  injection  time  corresponding  to  the  maximum 
photographic  area: 

+  sing),)  —  —  lg«tg<5«  )  ) 

J  («<90*-|<5®|)  (W) 

f»«*!,‘g~,OgSfr  3ia<P,')+  Mn“'(-tgttgd«)  j 
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or 


f„  =  12*  +  tg~'(  tg  0  sin  <p,)—  sin  -,(  —  ctg  i  ctg(3®  )  ) 

>  (i>90--j<5®i)  (9B) 

/«  =  tg_,(tg0  sinip,)  +  sin~‘(  —  ctgic:g<5®)  J 

2.  The  Situation  of  Observation  Satellites  With  Mission  of 
Photograph  in  Long  Flights 


At  present,  the  flight  times  of  observation  satellites  with 
missions  of  photograph  can  be  lengthened  to  several  weeks  to 
several  months.  This  type  of  satellite  generally  travels  on  a 
sun- synchronous  orbit.  By  appropriately  selecting  the  launch 
date  and  launch  time,  the  satellite's  photographic  area  can 
approach  global  coverage. 


As  regards  the  sun-synchronous  orbit,  - Q in  formula  (3B) 
can  be  considered  to  be  constant.  It  is  predetermined  by  the 
satellite's  flight  life,  solar  light  irradiation  conditions  of 
the  photographic  area,  the  satellite's  temperature  control  and 
other  requirements.  By  letting  gUQ-  Q=r,  seeking  the  deriva¬ 
tive  of  tg  for  6©  and  letting  its  derivative  be  zero,  we  can 
obtain: 

tg<$0  =  -  ( ioa) 

The  sun's  right  ascension  and  declination  &$>  have  the  fol¬ 


lowing  relational  formula: 

sinfl,0  = 


tg  t$6> 
tg  £ 


(10B) 


In  the  formula,  €  is  the  included  angle  of  the  ecliptic  surface 
and  the  earth's  equatorial  surface.  Therefore,  it  corresponds 
to  obtaining  the  launch  date  of  the  maximum  photographic  lati¬ 
tude.  We  can  find  £2  of  the  corresponding  injection  time  from 
the  following  formula  and  after  substituting  it  into  formula 


(5B)  we  can  obtain: 


(11) 


By  substituting  formula  (10A)  into  formula  (3B) ,  we  can 
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obtain  maximum  photographic  latitude 


max 


tg  =+tgi 

*  max  —  ’ 

Because  sun-synchronous  orbit  i>90°,  the  corresponding  photo¬ 
graphic  latitude  range  is : [180 °-i, i-180 0 ] . 

VII.  The  Relationship  Between  Solar  Altitude  and  Local  Time 

When  satellites  with  mission  of  visible  photograph  take 
pictures  of  earth  surface  targets,  it  is  sometimes  necessary  to 
take  the  pictures  at  specific  solar  altitudes  (the  altitude  of 
the  sun  from  the  horizon  is  the  elevation  angle).  Therefore,  it 
is  necessary  to  know  the  relationship  between  the  solar  altitude 
and  local  time . 

If  the  latitude  of  target  area  M  is  the  altitude  of  sun 

S  is  6  as  shown  in  Fig.  4. 


Fig.  4 

By  passing  the  M  point  and  making  a  declination  circle,  we  can 
obtain  spherical  surface  triangle  MOS.  From  this  spherical  sur¬ 
face  triangle  we  can  obtain: 

sin#  =sin  <S@  sin<PM-cos  cos^j^cos^  (12) 

The  above  formula  gives  the  relationship  between  solar  altitude 
0  and  local  time  t... 
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VIII.  Numerical  Examples 


Example  1.  A  certain  known  satellite's  orbit  insertion  point 
latitude  <T>z=35°  and  its  orbit  entering  direction  'f' = 37°.  We 
attempt  to  determine  the  injection  time  and  cause  the  ground  area 
which  can  be  photographed  by  the  satellite  to  be  maximum. 

Based  on  the  orbit  insertion  point  and  orbit  insertion  dir¬ 
ection  of  the  satellite,  we  can  find  inclination  i  of  the  satel¬ 
lite's  orbit  surface  by  the  following  formula: 

cosi=sin  y  cos  S3 

z 

By  substituting  the  data  of  the  satellite's  orbit  insertion 
point  latitude  and  orbit  insertion  direction  into  the  above  form¬ 
ula,  we  can  obtain: 

i=60 . 4  0 

Therefore,  we  can  obtain  the  two  injection  times  from  formula  (9A) 

tzl*=13. 5h-sin-1  (-1.76tg(5ff>  ) 
tz2=1.5h  -sin  1  (-1. 76tg($fi>  ) 

In  the  above  formula,  t  .  and  t  «  change  with  the  launch  date  as 

Z  i  Z  A 

shown  in  Fig.  5. 

* 

<+**i  (1) 


Key:  (1)  Hours;  (2)  Months. 


We  can  see  from  Fig.  5  that  this  satellite  can  obtain  two 
injection  times  corresponding  to  the  maximum  photographic  area: 
one  is  daytime  and  the  other  is  night. 

Example  2.  It  is  necessary  that  the  satellite  mentioned  in 
the  above  example  photograph  a  certain  focal  target  and  the  lat¬ 
itude  of  this  target  be  56°  north.  It  is  also  necessary  that  the 
local  time  for  photographing  this  target  area  be  10  A.M.  and  that 
we  try  to  determine  the  injection  time  of  the  satellite. 

By  eliminating  AQ  we  can  obtain  the  injection  time  from 
formula  (6)  : 

t  =7.6h 
z 

which  is  7:36  A.M. 

IX.  Discussion 

1.  The  launch  time  of  the  satellite,  aside  from  the  illumina¬ 
tion  conditions  of  the  subsatellite  point,  is  usually  also  limited 
by  many  other  factors.  For  example,  the  requirements  of  the  alti¬ 
tude  control  system,  the  solar  cell  system  and  the  temperature 
control  system  etc.  Therefore,  the  launch  time  of  the  satellite 
always  requires  compromise  determination  after  integrating  the 
requirements  of  many  areas.  However,  this  type  of  compromise 
consideration  should  be  established  on  the  basis  of  the  detailed 
theoretical  analysis  of  each  area. 

2.  If  the  satellite  is  required  to  photograph  specific  areas 
at  a  given  local  time  or  at  a  given  solar  altitude,  at  this  time 
we  can  calculate  the  injection  time  of  the  satellite  (example  2) 
from  formula  (6).  However,  for  this  reason  the  calculated  in¬ 
jection  time  is  the  only  necessary  condition  for  photographing 
specific  areas  in  a  given  time  but  it  is  not  the  full  condition. 

In  order  to  cause  the  satellite  to  pass  just  above  the  specific 
area  at  a  given  time,  it  is  also  necessary  to  appropriately 
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regulate  the  orbital  period  of  the  satellite. 

We  would  like  to  express  special  thanks  to  comrade  Fan  Qinhong 
for  his  reading,  revising  and  valuable  opinions  of  this  article. 
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SOCIETY  DEVELOPMENTS 


The  Dynamic  Environment  Conditions  Symposium  Convenes  in  Dalian 

The  "Structural  Strength  and  Environmental  Engineering" 

Special  Committee  of  the  China  Aerospace  Society  convened  the 
"Dynamic  Environment  Conditions  Symposium"  in  Dalian  on  August  17, 
1981.  53  representatives  from  23  units  participated  in  the  con¬ 

ference  and  19  technical  reports  were  exchanged.  The  conference 
was  chaired  by  special  committee  chairman  Zhang  Jingrong. 

The  professional  participants  joined  together  to  summarize 
China's  experience  in  dynamic  environment  engineering  of  carriers 
and  satellites,  analyze  the  present  situation  and  find  certain 
existing  problems  influencing  the  development  of  environmental 
engineering  such  as:  1.  the  knowledge  of  vibration  sources  is  not 
uniform  nor  are  data  processing  and  the  analysis  methods;  2. 
there  is  little  effective  telemetering  data  and  ground  measurement 
data  which  has  caused  blindness  in  conditions  formulation  work; 

3.  the  systems  level  test  conditions  are  deficient;  4.  not 
enough  attention  has  been  given  to  fundamental  research.  After 
serious  discussions  and  exchanges,  they  clarified  certain  long¬ 
standing  vague  ideas  and  brought  forth  seven  proposals  for 
improving  environmental  conditions  work.  There  was  especially 
unified  understanding  of  the  recent  pressing  need  to  develop  work 
and  they  proposed:  1.  to  actively  prepare  to  formulate  dynamic 
environmental  standards  for  carriers  and  satellites;  2.  to  inte¬ 
grate  the  development  of  new  models  to  carry  out  research  of 
environmental  control;  3.  five  proposals  to  develop  reliable 
research  to  integrate  with  environmental  engineering  and  research 
on  the  basic  theories  of  environmental  engineering. 

The  environmental  conditions  work  on  spacecraft  directly 
influences  its  performance  and  flight  reliability.  The  partici¬ 
pants  at  the  conference  unanimously  considered  that  each  level  of 
leadership  should  give  full  attention  to  this  work  and  all  at  the 
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conference  passed  the  "certain  major  problems  and  proposals  in 
formulating  spacecraft  dynamic  environment  and  test  conditions.” 

The  conference  triumphantly  concluded  on  August  22  and  Li 
Shuou  delivered  the  closing  speech.  He  said  that  environmental 
conditions  work  was  started  from  scratch,  that  great  achievements 
were  gained  and  that  now  was  the  time  for  summarizing  experiences. 
He  hoped  that  on  the  basis  of  this  conference  comrades  would 
further  summarize  experiences  and  write  articles  in  each  field  of 
specialization .  He  especially  pointed  out  that  the  formula¬ 

tion  work  of  environmental  conditions  was  the  central  link  in 
environmental  work.  He  hoped  that  everyone -would  draw  on  col¬ 
lective  wisdom  and  absorb  all  useful  ideas  and  as  quickly  as 
possible  formulate  feasible  standards.  He  also  said  that  we  must 
strengthen  research  work  on  the  environment,  replenish  the  tech¬ 
nical  contingency,  increase  necessary  equipment,  handle  well 
preparatory  work  and  make  deserved  contributions  to  research  on 
new  products. 

The  Astronautic  New  Welding  Techniques  Exchange  Meeting  Convenes 
in  Qingdao 

The  Materials  Technology  Special  Committee  and  the  Spacecraft 
Manufacturing  Technology  Committee  of  the  China  Aerospace  Society 
jointly  sponsored  the  first  "Astronautic  New  Welding  Techniques 
Exchange  Meeting”  in  Qingdao  from  August  21  to  25  of  1981. 
Altogether  69  representatives  from  27  units  participated  in  the 
meeting. 

The  aim  of  this  meeting  was  to  summarize  and  exchange  exper¬ 
iences  for  the  development  of  models  by  welding  techniques, 
exchange  results,  have  mutual  study  and  promote  the  development 
of  astronautical  welding  techniques. 

The  meeting  issued  a  total  of  40  articles  and  technical 
reports.  The  major  contents  of  these  articles  and  technical 
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reports  included:  materials  welding  performance;  astronautical 
structural  parts  welding;  braze  welding;  new  technical  equipment 
for  welding;  special  welding  as  well  as  welding  quality  control 
etc.  For  example:  all  of  the  articles  showed  that  many  of  the 
units  had  done  a  great  deal  of  work  on  the  welding  of  new  mater¬ 
ials  such  as  high  strength  aluminum  alloys,  ultra-high-strength 
steel,  martensitic  rust-proof  steel,  heat  resisting  alloys, 
integrated  product  structures  to  resolve  the  many  technical 
problems  of  product  welding  as  well  as  braze  welding  techniques, 
new  technical  equipment  for  welding  and  welding  quality  control 
and  analysis.  There  was  also  a  large  amount  of  work  done  in  the 
area  of  special  welding  techniques.  The  participants  at  the 
meeting  considered  that  after  exchange,  there  were  gains,  raising 
of  standards  and  predetermined  goals  were  reached. 

The  meeting  also  invited  Tian  Yitang  and  Wu  Lin  of  Harbin 
Engineering  University,  He  Fangdian  of  Qinghua  University,  Liu 
Yisheng  of  Harbin  Welding  Institute  etc.  to  present  reports  on 
welding  joint  fracture  toughness,  the  use  of  microprocessors  in 
welding  as  well  as  automatic  control  in  the  welding  process. 

These  reports  will  have  very  good  guiding  effects  for  future 
work.  The  meeting  considered  that  in  order  to  raise  astronautical 
welding  quality  and  raise  the  reliability  of  welded  components,  we 
should  carry  out  further  work  in  the  areas  of  welding  cracks, 
welding  fracture  toughness  and  automatic  control  in  the  welding 
process. 

The  meeting  also  held  discussions  on  how  to  hold  specialized 
academic  meetings  as  well  as  how  to  chose  academic  papers.  They 
also  proposed  certain  beneficial  recommendations.  Moreover,  it 
was  decided  that  the  "Selection  of  Papers  on  New  Welding  Techni¬ 
ques  in  Astronautics"  would  be  published  in  the  future. 

The  Special  Committee  on  Measurements  and  Tests  of  the  China 
Aerospace  Society  is  Established 


58 


The  Special  Committee  on  Measurements  and  Tests  of  the  China 
Astronautical  Society  convened  an  inaugural  meeting  and  its  first 
academic  exchange  meeting  in  Beijing  from  September  3  to  8  of 
1981.  Participating  at  the  meeting  were  88  representatives  from 
40  units  of  the  National  Defense  Science  Committee  System,  the 
Chinese  Academy  of  Sciences,  the  Chinese  Scientific  Research 
Institute  of  Measurements,  institutions  of  higher  learning,  the 
various  ministries  of  machine  building  etc. 

This  meeting  received  a  total  of  42  papers  and  reports.  Com¬ 
mittee  Chairman  Qiao  Shijing  presented  a  report  entitled  "The 
Place  and  Effects  of  Measurements  and  Tests  in  the  Aerospace 
Industry."  At  the  meeting,  representatives  separately  presented 
reports  on  "The  Important  Effects  of  Measurements  and  Tests  in 
Space  Test  Control  Techniques,"  "Developments  in  Measurements," 
"The  Use  of  Lasers  in  Measurements,"  "The  Automation  of  Measure¬ 
ments  and  Tests,"  "Error  Analysis"  and  "Dynamic  Parameter 
Measurements  and  American  Dynamic  Test  Techniques."  These  reports 
caused  the  representatives  to  understand  the  present  situation  and 
development  techniques  of  domestic  and  foreign  measurements  and 
tests,  broadened  their  outlooks  and  clarified  their  directions. 

Aside  from  these,  the  five  specialized  groups  of  overall 
planning,  long  thermodynamics,  radio  and  time  frequency  separately 
exchanged  papers  and  reports  and  there  were*  rich  contents  related 
to  each  specialized  field  of  measurements  and  tests  of  the  aero¬ 
space  system.  The  main  contents  were  foreign  and  domestic  dynamic 
states,  scientific  research  results,  automatic  testing,  data 
processing  and  error  analysis  as  well  as  research  on  measurements 
and  test  equipment  and  precise  measuring  methods  to  develop 
urgently  needed  products.  Pleasing  results  were  attained  in 
studying  the  close  integration  of  measurements  and  tests  and 
astronautical  scientific  research  and  production  and  in  resolving 
key  problems  of  products  and  special  measurements  and  tests.  In 
order  to  guarantee  the  precision,  quality  and  reliability  of  the 

59 


aerospace  system  we  brought  the  important  technical  foundations 
and  obstructions  into  play. 

In  the  exchanges,  the  representatives  reviewed  the  exper¬ 
iences,  lessons  and  achievements  in  astronautical  measurement  and 
test  work  over  the  last  twenty-five  years  and  fully  determined 
that  measurement  and  test  work  contributed  to  China’s  aerospace 
industry.  In  the  statements  of  many  of  the  representatives,  it 
was  proposed  that  measurement  and  test  work  should  run 
through  the  entire  process  of  advance  research,  test  manufactur¬ 
ing,  tests,  production  and  use  of  products.  It  is  one  important 
key  in  ensuring  the  quality,  reliability  and  even  the  success 
and  failure  of  products.  Many  of  the  attending  representatives 
hoped  to  draw  the  attention  of  each  level  of  leadership  and  re¬ 
lated  departments  because  of  the  instances  wherein  measurement 
and  test  work  was  unable  to  keep  up  which  influenced  the  progress 
of  product  research.  The  representatives  said  with  deep  feeling: 
"Measurement  and  test  work  must  come  before  product  development." 
In  view  of  the  fact  that  future  demands  for  measurement  and  test 
work  will  be  even  higher,  the  representatives  made  several  recom¬ 
mendations:  1)  it  was  recommended  that  an  assistant  designer  be 

in  charge  of  quality  and  measurement  and  test  work:  2)  we  must 
bring  forth  requirements  and  formally  take  on  responsibility  in 
the  product  planning  and  demonstration  stages;  3)  we  urgently 
need  to  reverse  the  trend  of. not  checking  special  testing  equip¬ 
ment  for  a  long  time,  formulate  systems  and  develop  work;  4)  we 
must  strengthen  the  departments  of  measurements  and  especially 
build  up  the  basic  units  of  the  industry  as  well  as  try  hard  to 
develop  advance  research  work  on  measurement  and  test  techniques 
and  fundamental  theories. 

The  meeting  also  formulated  activity  plans  for  each  special 
group  for  the  1982-1984  period  and  proposed  reliminary  ideas  for 
each  group  to  lead  the  units. 

The  representatives  considered  that  this  meeting  was  the  first 
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grand  meeting  of  the  specialty  of  astronautical  measurements  and 
tests.  The  astronautical  measurement  and  test  techniques  possess 
the  special  features  of  using  a  combination  of  many  branches  of 
learning.  This  meeting  established  exchange  channels,  strength¬ 
ened  connections,  had  mutual  study  and  relatively  large  achieve¬ 
ments.  Because  the  main  object  of  this  meeting  was  inauguration, 
time  for  exchange  was  short  and  the  number  of  participants  was 
small.  Thus,  it  could  not  completely  reflect  the  level  of  this 
specialty.  It  was  proposed  at  the  meeting  that  in  the  future 
special  groups  should  be  central  in  carrying  out  activities  and 
there  should  be  full  participation  by  those  in  the  same  field. 

As  regards  the  future  work  of  measurements  and  tests  and  the 
special  committee,  the  representatives  at  the  meeting  stressed 
that  in  order  to  adapt  to  the  demands  for  new  products  with  high 
precision  and  technical  indices,  measurements  and  tests  must  use 
various  new  techniques  and  develop  in  the  directions  of  automa¬ 
tion,  wide  frequency  bands,  wide  ranges,  many  functions,  dynamic 
tests  etc.;  strengthen  basic  theoretical  research  work;  strengthen 
the  training  of  measurement  personnel;  strengthen  the  integration 
of  measurement  and  testing  with  the  products;  handle  well  the 
journal  "Measurement  and  Test  Techniques  in  Astronautics"  and 
strive  to  make  larger  contributions  to  the  development  of  China’s 
measurement  and  testing  industry. 

Standards  for  Soliciting  Contributions  for  the  "Journal  of 
Astronautics 

A.  The  Journal  of  Astronautics  is  a  comprehensive  scientific 
journal  directed  by  the  China  Aerospace  Society.  It  uses  Marxism 
and  the  thought  of  Mao  Zedong  as  the  guiding  principles,  reflects 
China's  achievements  in  aerospace  science  and  technology,  promotes 
domestic  and  foreign  scientific  exchange  and  undertakes  to  realize 
the  four  modernizations. 

B.  The  Journal  of  Astronautics  is  a  quarterly  journal 


distributed  domestically  and  abroad.  The  readers  are  mainly 
those  engaged  in  scientific  research  in  the  aerospace  industry, 
engineers  and  technicians,  teachers  at  universities,  special¬ 
ized  schools  and  institutes  and  high  level  students  at  those 
institutions . 

C.  The  contents  of  the  journal  are  mainly  concerned  with  the 
theories  of  space  flight,  astronavigational  engineering,  techni¬ 
ques  (such  as  spacecraft,  power  devices,  control  systems,  launch 
installations  testing  instruments  etc.),  space  physics,  aero¬ 
space  medicine,  space  biology  and  space  law.  It  includes: 

1.  Certain  academic  reports  in  theory,  practice  or  the  two 
combined. 

2.  Technical  reports  concerning  achievements  and  advancements 
in  development  and  practice  as  well  as  those  which  reflect  pre¬ 
sent  advanced  levels. 

3.  Comprehensive  reviews  of  certain  theories  or  techniques. 

4.  Reports  and  comments  on  domestic  and  foreign  aerospace 
science  trends. 

D.  Requirements  for  Contributed  Articles 

1.  The  thesis  must  be  clear,  the  writing  concise,  the  data 
reliable  and  the  charts  and  tables  clear;  scientific  articles 
are  limited  to  10,000  characters  and  technical  reports  must  not 
exceed  5,000  characters;  the  beginning  of  the  article  must  have 
an  abstract  of  the  contents  (limited  to  300  characters)  and  an 
English  translation  of  it,  the  author's  name  must  be  written  in 
Chinese  pinyin  romanization  and  the  translation  and  Chinese 
pinyin  romanization  must  be  written  in  block  letters. 

2.  There  should  be  one  original  and  two  copies  of  the  contrib¬ 
uted  article  and  the  final  copy  should  be  a  clean  copy  according 
to  the  requirements  of  public  publishing;  it  should  be  written 


with  a  pen  on  one  side  of  squared  paper  and  each  character  and 
punctuation  mark  must  occupy  one  square;  the  writing  must  be 
clear,  script  characters  are  not  allowed  and  simplified  charact¬ 
ers  are  to  be  used. 

3.  Please  write  the  letters  inside  and  outside  the  manu¬ 
script  in  block  letters:  where  italics  are  necessary,  these  can 
be  written  in  script.  When  there  are  different  types  of  writing, 
please  write  notes  clearly  on  the  side  in  pencil  for  alphabetic 
letters  which  are  easily  confused  or  those  with  sizes  which  are 
not  easily  distinguished;  the  position  and  height  of  upper  and 
lower  corner  letters  and  numbers  on  arithmetic  symbols  must  be 
clearly  distinguished. 

4 .  The  charts  inserted  in  the  article  should  accord  with 
charting  standards,  use  prepared  Chinese  ink  to  draw  oh  transpar¬ 
ent  or  charting  paper  and  it  is  best  for  the  size  of  the  paper  to 
be  double  that  used  by  the  journal;  the  writing  in  the  charts 
should  be  done  in  pencil,  the  position  of  the  chart  in  the  manu¬ 
script  should  be  marked  by  a  square  and  the  number  of  the  chart 
and  its  title  should  be  given  below  the  square;  photographs  of 
objects  must  be  sharply  demarcated  in  black  and  white  and  be  clear; 
the  charts  and  photographs  should  be  separately  appended  and  not 
pasted  on  the  square  paper. 

5.  The  references  appended  at  the  end  of  the  paper  should  be 
from  well  known  domestic  and  foreign  books,  journals  and  articles 
and  it  is  requested  that  unpublished  material  not  be  used.  The 
writing  form  is:  number,  author's  name,  name  of  journal,  the  date, 
volume  and  number,  the  page  number  and  the  year. 

6.  It  is  requested  that  the  real  name,  work  unit,  address  and 
telephone  be  given:  please  do  not  submit  a  manuscript  many  times 
and  do  not  submit  articles  which  have  already  been  published  (with 
the  exception  of  restricted  publications)  for  reprinting;  for 


articles  published,  remunerations  will  be  given  according  to 
national  stipulations;  the  editorial  department  is  responsible 
for  returning  manuscripts  which  have  not  been  published. 

7.  Send  contributed  articles  to:  P.O.Box  838,  Editorial 
Department  of  the  China  Aerospace  Society,  Peking. 

Notice  of  Address  Change  of  the  China  Aerospace  Society 

Starting  January  1,  1982,  the  China  Aerospace  Society  will  have 
its  office  at  a  new  address. 

Address:  No.  2  Beixiao  St., 

Juwai  Yuetan,  Beijing 

Communications:  P.O.Box  838,  Beijing 

Telephone:  89-4602 

89-0651  ext.  303  (office) 

ext.  450  (editorial  office) 

Cable  address:  CSAPK 

The  original  P.O.Box  848  and • telephone  89-5040  will  be  specially 
used  by  the  China  Aerospace  Society  for  outside  liaison  depart¬ 
ments. 
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